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The Nambu Jona-Lasinio Lagrangian offers an explication of the seemingly contradictory obser- 
vations that a) the energy loss in the entrance channel of heavy ion reactions is not sufficient to 
thermalize the system and that b) the observed hadron cross sections are in almost perfect agreement 
with hydrodynamical calculations. According to this scenario, a critical opacity develops close to 
the chiral phase transition which equilibrates and hadronizes the expanding system very effectively. 
It creates as well radial flow and, if the system is not isotropic, finite V2 values. 

I. INTRODUCTION 

Some of the most surprising results of the ongoing RHIC experiments are the large radial flow 
Prad ~ 0.5, the large v-i values, where vi is the second Fourier coefficient of a Fourier expansion 
in (j) p of dN I (dyptdp t d(j)p) , and the fact that hydrodynamical models describe not only the total 
but even the differential anisotropy V2(pt) up to a transverse momentum of pt — 2 GeV [1-5]. The 
initial energy densities, at r = 1/m/c, employed in these calculation are between 4 and 6 GeV/ fm 3 
and agree quite well with the Bjorkcn estimate [6]. This seems to indicate that the system is at 
that time close to local equilibrium. For a survey of the results of the hydrodynamical approaches 
see [7]. The problem with these findings is that one does not understand how the system comes 
so fast to equilibrium. The energy loss of incoming partons is of the order of a 1 GeV/fm [8, 9]. 
Independent of the reaction scenario employed (minijets[10] or Weizsaecker- William fields of the 
incoming nuclei which materialize into on-shell gluons [11]) the time needed for equilibration is 
much longer. For a recent review see [12]. 

All these models share the feature that at energies below the validity of pQCD phcnomcnological 
approaches have to be employed. Therefore it may be worthwhile to use another phcnomcnological 
model which has been shown to be quite successful in regions where it can be compared with QCD, 
the Nambu Jona Lasinio (NJL) model[13, 14]. Supplemented by a t'Hooft determinant the NJL 
Lagrangian has the same symmetries as the QCD Lagrangian which are known to be the basis 
of many properties of quark and hadronic matter. It predicts[15], as pQCD [16, 17] calculations, 
color flavour locking at low temperatures and high densities, allows to describe the meson masses 
at low density and low temperature [14] and provides a simple approach to study the chiral phase 
transition. It predicts as well the tricritical point in the p — T phase diagram which has recently 
been predicted by lattice QCD calculations[18]. Of course this model has as well its deficiencies: 
Based on a local 4 - Fermion interaction it is not renormalizable and therefore a cutoff A has to be 
employed to regularize the loops. Furthermore due to the locality of the interaction confinement 
is not present and gluons do not appear as degrees of freedom. Due to these drawbacks it is hard 
to judge the quantitative prediction of this model but it can certainly serve for qualitative studies. 

Using this model we will show that due to the interplay between the quark and meson masses 
close to the chiral phase transition the system may develop a critical opacity where all s - channel 
transition rates become very large. This it true for the elastic as well as for the hadronization cross 
sections. These large cross sections equilibrate very efficiently the expanding plasma and create 
a gaz of mesons although confinement is absent in the NJL Lagrangian. Since the equilibration 
takes place during the expansion a strong radial flow develops which is seen the data and which is 
as well not understood yet. 
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II. THE NJL APPROACH 



The NJL model is the simplest low energy approximation of QCD. It describes the interaction 
between two quark currents as a point- like exchange of a perturbative gluon [13]. This local 
interaction is given by 

N c 2 -1 3 
c=l i,j 

where we have explicitly shown the flavor (i,j) and color/Dirac (a, /3, 7, S) indices. We normalize 
Si=o KtB^Sa — 2- Applying a Fierz transformation in color space to this interaction the La- 
grangian separates into two pieces: an attractive color singlet interaction between a quark and 
an antiquark (£(q q )) and a repulsive color anti-triplet interaction between two quarks £( qq ) which 
disappears in the large iV c limit. Here we are only interested in the color singlet channel (the color 
octet channel gives diquarks and can be used to study baryons [15]): 

£ = Co + £(qq) + £a (2) 

where £0 is the free kinetic part. Concentrating on the dominant scalar and pseudo scalar part in 
Dirac space we find the following Lagrangian which we use here: 



/={u,d,s} 



a=0 



(3) 



-G D [det[q f {\ - 175)9/ ] + det[q f {\ + 175)?/]] • 



The first term is the free kinetic part, including the flavor dependent current quark masses 
mo/ which break explicitly the chiral symmetry of the Lagrangian. The second part is the 
scalar/pseudoscalar interaction in the mesonic channel. It is diagonal in color. We have added 
the six point interaction in the form of the t'Hooft determinant which breaks explicitly the Ua{1) 
symmetry of the Lagrangian. The det runs over the flavor degrees of freedom, consequently the 
flavors become connected. 

The model contains five parameters: the current mass of the light and strange quarks, the 
coupling constants Gd and Gs and the momentum cut-off A, are fixed by physical observables: 
the pion and kaon mass, the pion decay constant, the scalar quark condensate < qq > and the 
mass difference between 77 and 77'. We will employ the parameters set:m° = 4.75 MeV, m {) s = 
147 MeV, G s /A 2 = 1.922, G D /A 5 = 10, A = 708 MeV 

The temperature and density dependent masses of the quarks are obtained by reducing the above 
Lagrangian to a one particle Lagrangian by contracting the remaining field operators in all possible 
ways. The meson masses are obtained by solving the Bethe Salpeter equation in the qq channel. 
The details of both calculations are found in ref. [14]. 

With these parameter we obtain the meson and quark masses displayed in fig.l. 



III. CROSS SECTIONS 



If created in heavy ion collisions the quark gluon plasma will expand rapidly. Therefore, not the 
static properties of the theory but the cross sections between constituents become dominant. In 
the NJL model these cross sections can be calculated via a 1/N C expansion [19, 20]. 

We start out with the elastic cross sections. The qq — > qq cross section has no s-channel contri- 
bution and is of the order of some mb. It is not important for the phenomena discussed in this 
letter. The Feynman diagrams for the qq — > qq cross section is shown in fig. 2. The matrix elements 
in the t and s channel are given by 

-iM t = S CuC3 S C2 . C4 u(p 3 )Tu(p 1 )[Wf '(pi - p 3 )}v(p 4 )Tv(p 2 ) (4) 

+ S c u c 3 S C2 .c 4 u(P3)(H5T)u(p 1 )[W^(p 1 - p 3 )]v{p i ){i-i b T)v{p 2 ) 

-iM s (p2)T«(pi)[iX>f (p! + p 2 )\v{p A )Tu{p 3 ) (5) 

+ S Cuc2 S c ^ C4 v(p2)(ij5T)u(p 1 )\iV^(p 1 + p 2 )]v(p 4 )(ij 5 T)u(p 3 ) , 
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FIG. 1: Masses of the the pseudo scalar mesons and of the quarks in the NJL approach. 




FIG. 2: Feynman diagrams of q + q — > q + q. 



where Pi(p2) is the momentum of the incoming q(q) and ^3(^4) that from the outgoing q(q). The 
Ci are the color indices and T are the isospin projections on the mesons. T> s and T> p are the meson 
propagators of the form 

V s ' p = ^ (6) 

1 - 2G s ns/p w 

with n s / p being the polarization tensor in the scalar/pseudoscalar channel. This cross section is 
displayed in fig. 3. We observe, as expected, an elastic cross section of the order of several millibarn 
but close to the critical temperature (T c — 24QMeV) the cross section increases dramatically due 
to the resonance structure in the s-channel. 

This resonance in the s-channel dominates also the hadronization matrix elements close to the 
threshold whose Feynman diagrams are given in fig. 4 

with the matrix elements: 

-iM t = ftSc lC2 v{p2)ilhigiSF{pi ~ Pz)ilhi92u(pi) (7) 
-iM u = fu5 Cl c 2 v(p2)i^igiS F {pi-pi)i^ig2u(pi) (8) 
-iM s = v(p2)S ClC2 f s iV(p 1 +p 2 )T(p 1 +p2;p3)igiig2u(pi). (9) 
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FIG. 4: Generic form of Feynman diagrams for t and s channel of q + q — ♦ M + M. The u channel is 
obtained by exchanging the mesons of the t channel. 



We discuss here as an example the hadronization cross section for u + u — * tt + + ir (which has 
two s and one t channel) and the cross section for the inverse reaction which are displayed in fig. 
5. 

Around T c the transition amplitude diverges close to the threshold, therefore the cross sections 
in both directions have their maximum there. The difference comes from the different flux and 
phase space factors. 



IV. EXPANDING PLASMA 



These cross sections suggest the following reaction scenario: as soon as projectile and target 
nuclei overlap the density is that high that the incoming nucleons overlap and partons can travel in 
this matter. By pQCD cross sections, which are of the order of some mb, the partons scatter but 
these cross sections are not strong enough to equilibrate the system. While the system expands 
the local available center of mass energy and hence the local temperature lowers. The closer the 
system comes to the phase transition the larger becomes the elastic cross sections. Due to this 
critical opacity the system behaves close to the phase transition more like a liquid then like a 
plasma. In this phase the system approaches thermal equilibrium. The large cross sections during 
the expansion create in addition a large radial flow and - if the system is azimuthally not isotropic - 
also V2 values. When passing the critical temperature the large hadronization cross sections become 
effective because the created mesons do not decay anymore. The inverse cross section is small due 
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FIG. 5: Cross section for u + u 
temperatures. 



7r + + tv and for the inverse reaction as a function of ^fs for different 



to the kinematic conditions. The system cools locally and the maximum of the cross section moves 
to larger values of ^fs. This means that quark pairs with larger center of mass energies can now 
hadronize. During the expansion the local temperature passes from T > T c to 0. Therefore the 
cross sections are large for a large interval of \fs as seen in fig. 6. 

There remains the quantitative question whether the strength of these cross sections is sufficient 
to hadronize the system completely: In order to check this we perform calculations using the 
Quantum Molecular Dynamics approach [22] which has been successful used to describe heavy 
ion reactions in this energy domain. In this approach partons are presented as Gaussians whose 
centers in momentum and coordinate space follow Hamilton's equations. In addition they interact 
via cross sections calculated from the same Hamiltonian. Fig. 7 presents a model calculation for 30 
quarks and 30 antiquarks and 100 pions initialized with isotropic distributions in coordinate and 
momentum space for an energy density of 3GeV/ fm 3 which give a temperature slightly above T c . 
There pions are unstable and therefore the mesons can decay into the lighter quarks. When the 
system approaches the phase transition the quark condensate < qq > increases: the quarks become 
heavier and this energy is taken out from the relative motion. At the phase transition the mesons 
become stable and lighter than the quarks. They escape from the system in which the quarks 
remain. By emitting pions the remaining quark system cools down until finally only 1 quarks and 
1 antiquark are left which have not found a partner to created two pions. Thus, despite of the fact 
that confinement is not enforced by the Lagrangian, under the condition of an expanding plasma 
confinement is almost complete. 

In conclusion we have shown, that the thermalisation observed in the spectra in relativistic 
heavy ion reactions may be reconciled with the small pQCD cross section assuming that the system 
passes during its expansion a phase transition with critical opalescence as predicted by the NJL 
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FIG. 6: Size of the cross section for u + u — > tv + + tv in the yfs - temperature plane. 
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FIG. 7: Number of pions and quarks as a function FIG. 8: Root mean square radius of the different 
of time. species as a function of time. 



- Lagrangian. Both, elastic as well as inelastic cross section, become large over a large kinematic 
region which is covered by the expansing system and lead to a local thermalisation as soon as the 
system comes to the phase transition. Because this transition takes places in an already expanding 
system, radial flow is created as observed experimentally. The observed v-i values are an image of 
the asymmetry at the beginning of the expansion. 
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